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CURRENT LARGE DEVIATIONS FOR ASYMMETRIC EXCLUSION 
PROCESSES WITH OPEN BOUNDARIES 

T. BODINEAU AND B. DERRIDA 

Abstract. We study the large deviation functional of the current for the Weakly Asymmetric 
Simple Exclusion Process in contact with two reservoirs. We compare this functional in the large 
drift limit to the one of the Totally Asymmetric Simple Exclusion Process, in particular to the 
Jensen- Varadhan functional. Conjectures for generalizing the Jensen- Varadhan functional to open 
systems are also stated. 
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1. Introduction 

Recently, a theory has been developed to describe the current fluctuations of diffusive stochastic 
gas in their non-equilibrium steady state [U El 03 03 HH1 132] • The large deviation functional of the 
current can be obtained in terms of a difficult variational problem; namely one needs to find the 
most likely time dependent density profile conditionally to the value of the current. So far some 
concrete informations have been derived: for example, explicit expressions of the current cumulants 
have been derived in J8] and dynamical phase transitions have been predicted in 13 IS] - I n the 
present paper, we analyze this functional in the case of a one-dimensional diffusive lattice gas of 
length N, the simple exclusion process, with a weak asymmetry v/N (this drift v/N represents the 
effect of an external field acting on the particles from left to right). We will in particular consider 
the large drift limit v — ► oo to investigate the relation between the hydrodynamic large deviation 
functional of the Weakly Asymmetric Simple Exclusion Process (WASEP) and the one of the 
Totally Asymmetric Simple Exclusion Process (TASEP) which was computed by Jensen- Varadhan 

EH- 

The density large deviation functional associated to the stationary measure of the TASEP in 
contact with reservoirs was computed in |161 117j . This computation was extended in [20] to the 
WASEP and the TASEP functional was recovered in the large drift limit (at least for some range 
of the parameters). This convergence was not a priori obvious as the hydrodynamic scalings in the 
WASEP and the TASEP are of different nature: for a system of size N, the time scales like N 2 in 
the WASEP and like N in the TASEP. Our goal here is to show that a similar convergence (from 
the WASEP to the TASEP) holds for the current large deviations for an open system in contact 
with two reservoirs (at least in some range of the values of the current). A similar convergence was 
already observed in |§] for the ring geometry. 

In the present paper, we also discuss the crucial role played by the weak solutions of Burgers 
equation as in the Jensen- Varadhan theory [21] . It is well known that the weak solutions of Burgers 
equation are not unique [SHI- A small viscosity enables to regularize the equation and to select the 
relevant weak solutions: from the PDE point of view, the others are disregarded as non-physical. 
These apparently non-physical weak solutions appear nevertheless as the optimal way to produce 
some current fluctuations and we show that the regularization by a small viscosity remains valid 
in the large deviation regime. We stress that the correspondence between the WASEP and the 
TASEP is obtained here at the macroscopic level (for the large deviation functionals). We are in 
fact not able to make a precise derivation starting from the microscopic models. 

l 
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The paper is organized as follows. The WASEP in contact with reservoirs is introduced in Section 
01 In Section 03 the current deviations of the WASEP are computed and in Section 03 they are 
compared to the Jensen- Varadhan theory which is outlined in Section 0J The generalization of the 
Jensen- Varadhan functional to open systems is also discussed in Section 03 

In the appendix we also calculate the probability of maintaining an anti-shock in the viscous 
Burgers equation, and in the formal limit v = N, we recover the Jensen- Varadhan functional |24| . 

2. Weakly Asymmetric Simple Exclusion Process 

In this section, we briefly review some known results on the hydrodynamic behavior of the Weakly 
Asymmetric Simple Exclusion Process (WASEP). 

2.1. The microscopic model. We consider the WASEP on the chain D^r = {1, . . . , N} in contact 
with two reservoirs at the boundaries. A particle configuration is given by rj E {0, 1} N , with rn = 1 
if a particle is at site i and r/j = otherwise. The WASEP is a Markov process evolving according 
to the following dynamical rules. For a given v > and N large enough, each particle attempts to 
jump to the right at rate \ — ^ and to the left at rate 5 + 277- If the target site is occupied, the 
jump is disregarded. 

Particles are injected or removed only at the reservoirs, i.e. at sites 1 and N. At the left 
boundary, i.e. at site 1, particles are created with rate cf and annihilated with rate c^. Similarly 
at the right boundary, particles are created at site N with rate and annihilated with rate c^. 

For any N, the expectation with respect to the invariant measure associated to the previous 
dynamics will be denoted by (-)n- 

Under the action of the dynamics, a particle current flows through the system. Let Q(i, [s,s']) 
be the integrated current through the bond (i, i + during the time interval [s, s'], i.e. the number 
of jumps from i to i + 1 minus the number of jumps from i + 1 to i. The space/time integrated 
current up to time s will be defined as 

N-l 

Q S = £G(*.[M)- (2.1) 

1=1 

Hydrostatic. 

We first consider the stationary regime. The steady state density will be denoted by p and the 
steady state mean current by q. At the macroscopic level, the domain Djv is mapped into D = [0, 1] 
and the local density of the steady state, at the macroscopic position x, can be defined as 



VxG]0,l[, P {x) = 




where (-)jv is the invariant for the measure. At the boundaries, p is given by the rates associated 
to the reservoirs 

p(0)=Pa = ^^ and p(l)=p b = + ■ (2.3) 

The local density satisfies Fick's law 

Vx€]0,l[, \&xP-W x (o{p)) =0, (2.4) 

with 

a(p) = p(l - p) . (2.5) 
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The hydrostatic equation (|2.4|) has been derived for a broad class of stochastic models |211 
HE El GDI and in particular for the WASEP in |2D|. 

The mean current can also be deduced from (|2.4|) (see 

lim lim ( — \ = --V x p + va(p) . (2.6) 



N^OO S^OO \ S / AT 2 

According to Fick's law the current through each bond is of the order 1/N, where N is the length 
of the system. The previous scaling comes from the fact that Q s is defined as the sum of the 
contributions of the N bonds (|2.1|h 

The mean current can be explicitly computed from (|2.t)|) . As an example, let us examine the 
case p a > 2 > Pb which becomes the maximum current phase in the TASEP, that is in the limit of 
a strong asymmetry [TJ1 OH E] • Equation (|2.6|) implies that 

Pa ^ l-Va I 

dp = / dv — 2 , 

where we used the change of variables p — > 1/2 + v with v a = p a — 1/2, = p — 1/2. This leads 
to the following asymptotics for large v 

9 =4 + ^ + °(^- (2 - 7) 

Hydrodynamic. 

We consider now the evolution of the stochastic process and denote by r/(s) the configuration 
at time s. Let ip be a macroscopic density profile. The evolution of the particle system over the 
time interval [0, s] and starting with initial data close to ip will be denoted by the probability IP-^ s 
(which depends on v as well) . More precisely, P^r s is such that the initial data are randomly chosen 
wrt the Bernoulli product measure with mean ip{i/N) at site i, then the time evolution is given by 
the Markov dynamics. To study the hydrodynamic limit, we introduce the diffusive scaling 

x = jf> t = ^- (2-8) 

At the macroscopic level, the particle system is identified to the macroscopic density profile p 
defined for any (x, t)elx [0, T] by 

(x+e)N 

p(x,t) = 1xm lim V rn{N 2 t) . (2.9) 

i=(x-e)N 

If initially the discrete system is close (after rescaling) to the macroscopic profile p(x, 0) = <p{x) 
then for any t > the rescaled system remains close to the density profile solution of 

V(x,i)6Dx[0,T], d t p = -Ap-vV x (a(p)) with p(0, t) = p a , p(l, t) = p b , (2.10) 

where <p is the initial data and u{p) = p(l — p). This holds with a probability F^ TN 2 close to 1 

HUIiaHEllHEl- 

For any (x, t) € ID x [0, T], a local version of the space-time current can be defined as a space-time 
average on the microscopic boxes {(x — e)N, (x + e)^} x [tN 2 , (t + e)N 2 ] 

(x+e)N 

i(x,t) = limJim o -2]f2 E Q{h[tN 2 ,(t + e)N 2 }). (2.11) 

i=(x—e)N 
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The scaling in (|2.11j) can be understood as follows: j is the rescaled current while the microscopic 
current through each bond is j/N. Then summing over 2eN bonds during a time interval eN 2 
implies that one has to normalize by a factor 2e 2 N 2 . The conservation of the number of particles 
at the microscopic level becomes through (|2.9j) and (|2.11j) 

d t p = -V x j . 

2.2. Large deviations. The hydrodynamic equation (|2.10j) describes the typical behavior under 
the stochastic evolution. The probability of observing a different trajectory can also be computed. 
In fact, one can even estimate the probability of observing a given current of particles j(x,t) and 
its corresponding density profile defined by 

V(x,t) GBx [0,T], d tP = -V x j, (2.12) 

with some initial data p(x, 0) = (f(x) in D (Note that the current cannot be completely arbitrary: it 
should be such that the density (|2.12|) remains in [0, 1] and is equal to p a and pf, at the boundaries 
at any time t > 0.) 

The probability of the event {r] ~ (j,p)} that the rescaled particle system has a current j and 
remains close to the density profile p introduced in (J2.12|) is given asymptotically in N by 

F n,tn* (iv ~ U> P)}) « exp ( - iVJf Q)T] (j, p) + o{N)) , (2.13) 

and the large deviation functional is 

„ , s f T f U(x,t) - ua(p(x,t)) + ±V x p(x,t)) 2 

Zfol&P) = dt dx^^ fy-1 L2 a . (2.14) 

1,1 Jo Jo 2 a(p(x,t)) 

This can be interpreted as a local Gaussian fluctuation of the current with a variance cr(p(x, i)J 
depending on the local density. 

The hydrodynamic large deviation theory was originally introduced in |27U38ll2*5] to estimate the 
probability of events corresponding to atypical space/time density profiles. Recent developments 
13 El 1251 132j led to the more general expression Q2.14JI which enables one to control a deviation 
of the density as well as the associated current (see the appendix of [5j for a heuristic derivation). 
Note that the current j(x,t) appearing in (|2. 131 - 0. 14|) is an average over a long microscopic time 
of order iV 2 (see ([2.11)1 V Recent works [181 119j have also considered the correlations between the 
density and the instantaneous current. At present we do not see any direct connection between 
their distribution of the instantaneous current and our j(x,t). 

3. Current large deviations for a system with reservoirs 

In this Section, we investigate the current large deviations for the WASEP in contact with 2 
reservoirs at densities p a and pb- 

3.1. General estimates. We are interested in the asymptotic probability of observing a deviation 
of the integrated current over a very long time interval [0,TiV 2 ], i.e. in the event that 

^- = i£dtj\xj(x,t) = J (3.1) 

for a very large macroscopic time T (the correspondence between Q and j was stated in ()2.11j) ). 
According to the large deviation principle (|2.13|) , one has to minimize the functional (|2.14|) over the 
time dependent density profiles under the mean current constraint (|3.1)1 J30IH]- Suppose that the 
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optimal way to achieve this deviation is by maintaining the density profile close to some optimal 
time independent density profile. Then, the variational problem simplifies and one has for T large 

with 

p [Jo 2(J {p{x)) 

where the infimum is taken over all the density profiles p(x) with boundary conditions p a , pb 
prescribed by the reservoirs. For simplicity, we wrote p' = V x p. 

This stationary regime assumption was already made in ;8j when v = and led to an explicit 
prediction for all the cumulants of the integrated current Qt in diffusive lattice gases. It was 
however understood in Ej that the stationarity assumption may not be satisfied for some diffusive 
lattice gases. In this case, a space/time dependent current is more favorable than a constant current 
and G{J) does not give the correct order for observing the mean current J = ^ J t=Q J B j(x, t) (see 

piainj). 

The expression (|3.2jl can nevertheless be interpreted as the asymptotic cost for observing a 
constant current j(x,t) = J for all times t in [0, T] (see 1)2.11(1 ). 




Let us first show that the optimal profile p of the variational problem (|3.2j) satisfies the relation 

(\A ={J-Mp)f + 2Ko(p) 1 (3.3) 

where K is an integration constant determined by the boundary conditions p a ,Pb (this formula is 
an extension of equation (15) obtained in [S]). To see this, we write for a given density profile p 

where H(p) = — ■ The corresponding Euler equation is 

w / r v / , WAP) n 
H{p) -Mp) + ^(pr = - 



Multiplying by p' we obtain 



This completes the derivation of (|3.3|) . 



In order to obtain more precise results, we are going now to specify the current deviation J as 
well as the boundary conditions. 
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3.2. Maximal current phase. As an example, let us consider the case p a > 1/2 > pj,, which 
becomes in the large v limit the "maximal current phase". According to (|2.7jl . the mean current is 
of order We are interested in the large v asymptotics of the large deviation functional G{yq) 
(given in (|3.2|l ) for current deviations of the form 

J = vq. 



Case 1 (current reduction): q < 1/4. 

Let a be the density for which the current is equal to uq in a stationary regime, we write 
q = o~(a) = <r(l — a) < 1/4 (with a > 1/2). As we will see the optimal profile to achieve this 
current deviation concentrates close to a or 1 — a as v tends to infinity. 

Setting K = u 2 k, the equation (j3.3|) of the minimizers of the functional G(vq) can be rewritten 



(q-cj{p)f + 2ka{p). 



(3.5) 



Different behaviors occur depending on the value of q, i.e. on the location of a wrt the reservoir 
densities p a and pb- 

• Suppose that a G [pb,p a ] or I — a £ [pb,p a ]- As the optimal density profile goes smoothly from 
Pa to pb, it has to take the value a (or 1 — a) and (|3.5j) implies that k ^ 0. This implies that the 
optimal profile is decreasing. Integrating (j3.5|) leads to 



1 



2v = / dp (3.6) 
J Pb ^(q-a(p)Y+2ka{p) 

Thus as v tends to infinity, k must vanish (in fact from (|3,6|) . one can show that k vanishes 
exponentially fast wrt v). The optimal profile remains essentially equal to a or 1 — a. 
The expression (|3.5|) of the optimal profile combined with ()3.2|) leads to 

pa 1 r ( q - <j{p)) 2 + ka(p) 



As k 



G{vq) 
when u — > 
lim 



dp 



2a(p) 



^{q-o-{p)Y+2ka{p) 



(Q ~ <r{p)) 



X: 



(q - o-{ P )) 2 + ko-(p) 



k^o y /( q - a ( p) )2 + 2 ko-(p) 




if p € [1 — a, a] , 
otherwise . 



(3.7) 



(3.8) 



Let us define p + = min{a,p a } and p~ = max{l — a,pb}, then only the densities in [p~ , p + ] 
contribute to the asymptotics 



lim —G{uq) 

v— »oo v 



°-{p) 



log 



p' 



1 



log 



p 



1 



(3.9) 



where we used that a(p) = p(l — p). 



In particular, for [1 — a, a] C [pb,p a ], this leads to the first main result of this paper 



lim —G{vq) 

v— >oo v 



where 



g(l - a) 

11 

g(u) = u(l - u)log' 



5(a)) , 



(3.10) 



(3.11) 



1-u 

In this case the optimal density profile is essentially a piecewise constant function with three jumps 
over regions of width \jv: a boundary layer near the left reservoir where the density varies from 
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p a to a, another one from 1 — a to pb near the right reservoir and a jump from a to 1 — a. The 
sharp discontinuity between a and 1 — a will be interpreted in Section 15.11 as a convergence to 
an anti-shock and the asymptotic cost (|3,1U|) coincides with the Jensen- Varadhan functional which 
will be defined in 

• Suppose that [p , p a \ C [1 — a, a]. 

The optimal profile has to be close to a or 1 — a otherwise G(vq) would scale like v 2 (see 1)3.4(1 ). 
thus the optimal profile cannot be always decreasing. We first assume that it initially increases 
from p a to a local maximum j u depending on v. According to (|3.5|) . one has 



= {q-a{ lu )Y + 2ka{ lv ). (3.12) 

The density cannot increase beyond 7^ otherwise the RHS of (|3.5|) would be negative. Thus the 
profile p(x) increases from p a to 7^ in the interval [0, x v ] and then decreases to pb in the interval 
[x u , 1]. The analogous of (|3.6|) can be rewritten as 



2V Xy 



dp 



1 



Pa 



^( q -a(p)Y+2ka{p) 



and 



2v (l-x u ) 



dp 



pi, 



^{ q -a{p)Y + 2ka(p) 



For large the important feature is the singularity at j u which implies that j u has to converge to 
a and k to as v diverges. In this case, the profile concentrates close to a. Then Q3.4I) becomes 



{q-a(p)) 2 + ka(p) 
^/{q-a{p)f + 2ka{p) 



+ (q - a{p)) 



dp 



(>i, 



(q - <j{p)) 2 + ka(p) 
y /(q-a{p)) 2 + 2ka(p) 



(Q ~ (rip)) 



In the limit v — > 00, a computation similar to (|3.9j) applies: the contribution of the reservoir p a 
vanishes and the asymptotic cost is 



cr(p) 



a 



p b + a(l - a) log 



Pb 



log 



o 



(3.13) 



Similarly, if the profile first decreases from p a to 1 
asymptotic cost is 



{I-Pb) ~~°{l-a) 
— a and then increases to pb, then the 



'l-a 



v(p) 



Pa 



[1 — a) — a(l — a) ( log 



Pa 



(1 " Pa) 



log 



(1 



OL, 



a 



(3.14) 



Therefore, we can state our second main result. If one defines the boundary cost by 

(u v \ 

log jz r "log 7- r , (3.15) 

then the optimal cost of the current deviation is given by 

lim -G{yq) = mm {F rcs (p a , I ~ a),- F ms (p b , a)} . (3.16) 

where as before q = a(a) = q(1 — a). We remark that the cost is due to the boundary effects and 
that the contributions from the left and the right reservoirs decouple: the first term in (|3,16|) is a 
cost for a boundary layer near the left reservoir with a density jump from p a to 1 — a whereas the 
second term is the cost for a boundary layer near the right reservoir with a density jump from a to 
Pb- One can trigger a phase transition between the optimal profile close to a and the one close to 
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1 — a by varying p a and pb (the critical curve is p a = 1 — pb). This provides an additional example 
of a non-equilibrium phase transition associated to a current large deviation [01 E] . 

Case 2 (current increase): q > 1/4. 

Once again the optimal profile is determined by ()3.5j) . The optimal profile decreases from p a to 
Pb and equation (|3.6|) applies. The minimum of the derivative p' (see (|3.5|1 ) is reached at p = 1/2, 
so that ()3.6j) implies that for large v 

2 



A; 



(3.17) 



where e v vanishes as v tends to infinity. The functional (|3.7|) can be rewritten, using 



G{vq) 
This leads to 



-u z k + v 



dp 



pi, 



Mp) 



^( q -a{p)Y+2ko{p) -{q- a(p)) 



G{vq) ~ lu 1 



q ~~4 



and 



lim —G{vq) 

v— >oo v 



oo . 



(3.18) 



A deviation q > 1/4 leads to asymptotics of G{vq) in v 2 in contrast to a deviation q £ [0,1/4] 
where G{yq) scales like za 

3.3. Boundary effects. The current deviations select very specific density profiles. We are going 
now to estimate the large deviations associated to atypical density profiles instead of atypical 
currents. The main motivation is to compute the boundary contribution for a large drift v. 

We first fix p a > 1/2 and let pb = 7. We are interested in evaluating the probability of main- 
taining a density profile close to 7 over a very long time T, i.e. 

P N,TN* ( Vt G [°» T ]> J dx (P( X ' *) " < £ ) « eX P ( " TNQ{l>, Pa , 7 )) , 

where the initial data is the profile uniformly equal to 7. The choice pb = 7 is made in order to 
suppress the boundary effect of this reservoir so that the main contribution will depend only on p a . 
As the constraint is time independent, the large deviations (|2.13j) imply that 

(j- vcr [ P {x))+y{x)f\ 



under the constraint that Jq dx (p( 
We are going to show that 




2a{p{x)) 



(3.19) 



0-7) 2 ^e- 



2' 



lim lim -G{v,p a ,l) 
e^O v— *oo v 



9(l) -^(l-T), 

0, 



if 7 e [0, 1 - p a ] 

if 7 G [1 - Pa, 1/2] 

if 7 G [1/2,1] 



(3.20) 



where F res was introduced in (|3.15|) and g in ()3.10() . This can be intuitively understood in view 
of the previous results on the current large deviations. F ies (p a ,~f) was computed in (|3.15|) as the 
optimal cost for producing a current vcr^) with an optimal profile close to 7 < 1 — p a , thus the 
constraint {f Q dx (p(x) — 7) 2 ^ e} does not change the asymptotic. If 7 € [1 — p a ,Pa] then the 
optimal profile to achieve a current vci^j) is concentrated close to 1 — 7 and 7 (see (|3.10jl ). In this 
case, the cost is essentially given by the sharp jump between 1 — 7 and 7. The location of the 
discontinuity is not relevant to compute the asymptotic cost 5(7) — g(l — 7)- To satisfy the density 
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constraint, it is enough to shift the discontinuity close to the reservoir p a so that all the energy 
becomes concentrated in the boundary layer. 

We turn now to the derivation of (|3.2U|) . The optimal current J (in the variational problem 
Ij3.19jl ) is related to the optimal profile p by 

j f 1 dx - / i j= ^w-w _ rl (3.2i) 

Jo cr{p(x)) Jo °(p(x)) J-y 2 a {p) 

The density constraint implies that for v large, the current is of the order J = vq, where q G 
1$ = Wil) ~ ^ a (l) + an< i ^ vanishes as e tends to 0. Therefore, one can substitute the density 
constraint by a current constraint and derive a lower bound for G{y,p a ,"f) 

S(,,A.,7)> inf infl/'fc ^-^ + Mx)) 2 ! 

l^h P [Jo 2a{p(x)) J <? e/ * 

where G was introduced in (\3.2\i . For 7 G [0, 1/2], then the asymptotics of G{vq) follow from Q3.1U[) . 
(|3.14j) . When 5 tends to 0, we deduce by continuity that the RHS of (j3.2Uf) is a lower bound for 
the asymptotics of Q(v, p a , 7). For 7 G [1/2, 1], a cancellation similar to the one in (J3.9|) occurs and 
we recover as a lower bound. 

In order to derive an upper bound for Q(y, p a ,7), it is enough to compute the asymptotic cost 
of well chosen density profiles. Let us simply consider the case 7 G [1 — p a ,Pa] and define 

1-A 2 = (a( 7 ) - <y(p)f + e'(p - 7 ) 2 + 2ka(p) , (3.22) 

where k is chosen such that the boundary conditions are satisfied. By tuning e' appropriately wrt 
e, we see that this profile concentrates to the density 7 for large v, so that the density constraint 
is satisfied 

2 a[p(x)) 

Proceeding as in (|3.7j) . H3.9|) . we deduce that the asymptotic cost of the profile ()3.22|) converges to 
5(7) — 5(1 — 7)- Similar arguments for the other values of 7 enables to complete the proof of (|3.2Uj) . 



If p a < 1/2, then similar computations lead to 



lim lim -Q(y,p a ,i) = < 



^ reS (Pa,7), 
F res (p a ,l- 7 ), 

0, 



if 7 G [0, 1/2] 

if 7 € [1/2,1 -p a ] , 

if 7 G [1 - p a , 1] 



(3.23) 



where F Tes was introduced in ()3.15|) . 

The boundary cost of the right reservoir can be deduced by symmetries (particle-hole, left-right). 
For p > 1/2, the asymptotic cost of a discontinuity 7 is given by 



lim lim —Q(f,j, pi,) = < 

e— >0 u^oo 11 



-F Ies (p b ,-y), 
-F res U,l-7), 



0. 



if 7 G [1/2,1] 

if 7 € [l-p 6 ,l/2] 

if 7 G [0,1 -pb] 



(3.24) 
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For p b < 1/2 



lim lim -Q(u,j,p b ) 

e— >0 v-^oo v 



-F^(p b , 7 ), 

5(1 - 7) - a{i), 
0, 



if 7 G [1 - Pb , 1] 
if 7 G [1/2, 1 - p 6 ] 
if 7 G [0,1/2] 



(3.25) 



In section l3~2l we focused on the current large deviation for p a > 1/2 > p b . For other choices of 
the reservoir densities, similar computation would lead to asymptotic costs which can be expressed 
also in terms of (j3.20j) . (j3.23M3~25j) . Further interpretation of the boundary effects is postponed 
until Subsection 15,21 

4. Jensen- Varadhan functional 

In this Section, we recall some properties of the TASEP and discuss the hydrodynamic large 
deviations which were derived by Jensen and Varadhan [21] . The connection between this functional 
and the previous results obtained for the WASEP will be detailed in Section El 

4.1. TASEP. Following the notation introduced for the WASEP, we consider the TASEP on the 
domain B^v = {1, • • • ,N} in contact with two reservoirs at the boundaries and denote a particle 
configuration at time s by rj(s) G {0, 1}^. For the TASEP, the particles only jump to the right 
with a constant rate 1 if the neighboring site on their right is empty. The reservoirs follow the 
same dynamical rules as the ones introduced in Section 12.11 Let be the probability measure 
associated to the TASEP on B^r during the microscopic time interval [0, s] and starting from the 
Bernoulli product measure with local density at site i equal to <p(i/N). 

The hydrodynamic behavior obeys a different scaling from the one of the WASEP 
At the macroscopic level, the particle system is identified to the macroscopic density profile 

(x+e)N 

p(x,t) = lim lim — — > riANt) . 

i=(x-e)N 

If initially the discrete system is close (after rescaling) to the macroscopic profile p{x, 0) then for 
any t > 0, it was proven in (see also [HHIEB]), that the rescaled system remains close to the 
density profile evolving according to Burgers equation 

V(s,t) GBx [0,T], d t p + V x (a(p)) =0 with a(p) = p(l - p) . (4.2) 

For general initial conditions, the solution of this PDE is not well defined as shocks may occur after 
a finite time in the bulk and at the boundaries. For simplicity, we will not discuss the boundary 
effects and focus on the bulk singularities. A detailed mathematical study of the PDE (|4.'2[) can be 
found in |2l 1361 15*7] . Notice also that the hydrodynamic limit can be derived beyond the shocks for 
a broad class of lattice gases ^ ESI • 

In order to take into account the shocks in the bulk, Burgers equation has to be rephrased in 
terms of weak solutions. Let $ be a smooth function with compact support strictly included in 
B x [0, T], i.e. there is 5 > such that 

[6,1-6], §(x,t)=0 and Vx G B, §(x,Q) = *(z,T) = 0. 
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Integrating 1)4. 2|) by parts leads to 

f dt I dx{pd t ^ + a(p) V x $} =0. (4.3) 
Jo Jo 

A trajectory p is said to be a weak solution of 1)4.2(1 if it satisfies (|4.3|) for any test function <I> with 
compact support strictly included in D x [0, T\. This formulation is not sufficient to prescribe the 
behavior at the boundary, nevertheless we will be mainly interested in bulk properties for which it 
is enough to consider (|4.3|) . 

For a given initial data, there is no uniqueness of the weak solutions. Let us consider, as an 
example, a weak solution on R. We define a density profile with a discontinuity a, b moving at 
speed v , i.e. 

fa, if x < vt , 

POM) = \ , - f ' (4.4) 

v ' [ 6, if x > vt . y ' 

This profile is a solution of 1)4. M|) for any a, 6 in [0, 1] provided that the velocity is given by 

v = a{a) ~ . (4.5) 
a — b 

If a < b, the previous density profile ()4.4j) is a physical solution. For a > b, the initial data 
with the discontinuity a, b at evolves as a rarefaction fan so that the previous density profile 
should be disregarded and considered as a non physical event. In the following, we will refer to the 
discontinuity a < b as a shock and a > b as an anti-shock |36| . 

A way to select the physical solutions is to add a small viscosity e and to recover Burgers equation 
in the limit e — > 

V(x,t) GDx [0,T], dtp + V^p)) = |Ap. (4.6) 

The previous equation can be identified to the WASEP evolution ()2.1Uj) after rescaling the time by 
a factor v = 1/e. Thus in the hydrodynamic regime, the TASEP can be seen as the limit of the 
WASEP at large drift. Formally, this means that at the hydrodynamic level the limits N — » oo 
and v — > oo can be replaced by v = N — > oo. 



4.2. Large deviations. The Jensen- Varadhan functional J\o,t]i introduced in |24| . provides an 
estimate of the probability that the rescaled particle configuration remains close to a given density 
profile p(x, t) over the macroscopic time interval [0, T] 

%,tn ({V ^ P}) * exp ( - NJ %T] (p)) , (4.7) 

with initial data </?(a?) = p(x, 0). The key role of the weak solutions ()4.3I) in the large deviation 
regime was understood in [21] : if p is not a weak solution then v7[o,r](/°) = oo, which means that the 
probability (|4.7|) vanishes at an exponential order faster than N. For general weak solutions, the 
Jensen- Varadhan functional [21] is defined in terms of the Kruzkhov entropy jSH]. In this paper, 
we will simply focus on a few concrete examples for which an explicit representation of i7[o,t] can 
be stated. 

Let p be a weak solution made of several shocks and anti-shocks (see 1)4.4)1 ) and such that p 
remains constant equal to p a and pi, in a small neighborhood of the reservoirs during the time 
interval [0, T] (the latter condition enables us to disregard the boundary effects). If p has a single 
anti-shock of the type ()4.4)) (i.e. a > b) during time T, then the Jensen- Varadhan functional derived 
in [21] is given by 

J [0tT] (p)=TF(a,b) 
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with 

F(a, b) = g(b) - g(a) - ^ ~ ^ (h(b) - h(a)) , (4.8) 

where 

u 

VuG[0, 1], h(u) = ulogu + (1 — u) log(l — u), g(u) = u(l — u) log u. (4.9) 

1 — u 

(Remark that h is the natural entropy associated to the TASEP and that g' = h'a' '.) If there are 
several anti-shocks their contributions add up: suppose that the evolution p contains t anti-shocks 
{(a,i,bi)}i <^e each of them maintained during the time interval C [0, T], then 

e 

J[o,T](p) = ^2(8i-ti)F(ai,bi). (4.10) 
i=i 



In the previous examples, the boundary effects are disregarded. In fact, the Jensen- Varadhan 
theorem was derived on a torus and does not take into account the influence of the boundaries. 
Nevertheless Q4.7JI should also apply in an open system for weak solutions p which are smooth 
at the boundary. Heuristically, the Jensen- Varadhan functional is given by the local cost of each 
anti-shock thus the regularity of the density at the boundary should be enough to decouple the 
bulk from the reservoirs. The boundary contribution will be discussed in Section [5.21 



5. From WASEP to TASEP 



5.1. Physical relevance of the weak solutions. Some of the weak solutions of the Burgers 
equation (eg. the anti-shocks) are considered as unphysical and disregarded. On the other hand, 
we showed in (|3.1(J|) that the anti-shocks with discontinuities (a, 1 — a) arise naturally as the 
limit of the optimal density profiles which realize a current large deviation va{a) when v — > oo. 
Furthermore, the large deviation functional associated to the WASEP converges to the Jensen- 
Varadhan functional corresponding to the stationary anti-shock with discontinuity (a, 1 — a) (see 
(j3.1U|) and l|4.8|) ). There are many ways of producing a current deviation even in a stationary 
regime, but conditionally to observing a current deviation of the form va(a), the typical profiles 
are close to an anti-shock with discontinuity (a, 1 — a): the current deviation selects a specific weak 
solution. 

Our observation shows that the viscosity regularization (|4.6|) remains valid beyond the law of 
large number regime and that the Jensen- Varadhan functional can be understood as a limit of the 
WASEP functional (|2.1U|) . Another example of an anti-shock traveling at non-zero velocity will be 
discussed in the Appendix. 

The correspondence between the WASEP and the TASEP process in the large deviation regime 
has been established at the macroscopic scale, i.e. for the large deviation functionals. At the 
microscopic level, the correspondence is far from being obvious as the scaling limit for the WASEP 
(I2.8j) and TASEP (|4.1|) are of different nature. We can rewrite (jlH.lOjl as 



uq 



exp 



-t 



>M + o(i 



exp 



tu 

"N 



F(a, 1 — a) 



(5.1) 



The order of the limits is the following: first t — > oo, then iV — > oo and finally v — > oo. At least 
formally, in the limit v = N — > oo, the asymptotics (|5.1I) is consistent with the Jensen- Varadhan 
theory. 
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Notice also that for the TASEP, the current distribution ^3 scales differently depending 
whether the current is decreased or increased. In (|3.18j) . a similar behavior was pointed out for the 
WASEP. For v large, the probability of observing a current deviation vq with q > 1/4 is given by 



N.t 



Qt 



vq 



cxp 



-t 



cxp 



'~N 




(5.2) 



At least formally, if v = N the previous asymptotics is consistent with the scaling expected for 
the TASEP 15 : to increase the current beyond 1/4 all the particles have to be accelerated. One 
should not however expect the leading order in (|5.2|) with v = N to be also valid for the TASEP. 
(|5.2|) implicitly assumes that for large N, the measure becomes locally a product measure. We 
believe instead that in the TASEP, the optimal measure to increase the current remains different 
from a product measure even for large N. 

5.2. Boundary contributions. The Jensen- Varadhan functional introduced in Section f4. 21 takes 
into account only the bulk large deviations and not the boundary contribution. At present, there 
is no microscopic derivation of the large deviations for the TASEP in contact with reservoirs. For 
diffusive models, the large deviation functional of the system with reservoirs is simply the functional 
in the bulk with the constraint that the densities at the reservoirs are fixed equal to p a and pb (see 
(|2.14l) or [3] for a rigorous derivation in the case of the SSEP). The Jensen- Varadhan functional 
cannot be generalized in this way, as we know that even in the stationary regime, the reservoirs do 
not impose a fixed density: the steady state of the TASEP can be discontinuous at the boundary 

O [23 ED El EH]- 

In Section 13.31 we computed the probability of observing a sharp variation of the density at the 
boundary for large v. As the Jensen- Varadhan functional can be interpreted as the limit of the 
WASEP large deviation functional in the bulk, we conjecture that (|3.2U|) . (|3.23M3~2*5|) provide also 
the right expression for the boundary contribution. If this is the case, the functional associated to 
the system with reservoirs should be the sum of the Jensen- Varadhan functional (in the bulk) and 
of the boundary effects. 



The previous conjecture cannot be validated without a proper derivation from the microscopic 
model. Nevertheless, further comments can be made to justify our guess. We will focus on the 
case p a > 1/2. For the TASEP on the half-line {1, 2, . . . }, one can find, for any density 7 in 1], 
an invariant measure with mean asymptotically close to 7 [HJ- This implies that there is no cost 
to maintain a discontinuity (p a ,j) for any 7 in [i, 1] as predicted by (|3.2l)j) . For 7 G [1 — p a , I], 
(|3.2U|) asserts that 1 — 7 plays the role of an effective boundary. The boundary layer is made of two 
pieces, first a jump from p a to 1 — 7 (which costs nothing as 1 — 7 > 1/2) and then an anti-shock 
with zero velocity between 1 — 7 and 7 which has a cost given by the Jensen- Varadhan functional. 
A similar interpretation holds for the expressions (|3.23M3~2*5|) . 



6. Conclusion 

In the present paper, we have analyzed the large deviation functional of the current for the 
WASEP in contact with reservoirs. In the large drift limit, we have obtained an expression (|3.1U|) 
similar to the Jensen- Varadhan functional (|4.8|) . We have also obtained the cost (13.20(1 , (|3.23M3~23|) 
of maintaining boundary layers in the large drift limit and we believe that these expressions should 
remain valid for the TASEP. Of course, it would be necessary to validate our expressions from a 
calculation starting from the microscopic model. For general weak solutions, the Jensen- Varadhan 
functional has to be formulated in terms of the Kruzkhov entropy [21]. It would be interesting to 
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relate the boundary costs (|3.2U|) . Q3. 23^3^21))) to the more general entropy condition introduced for 
open systems in [2"| B7|. 

A strong connection was derived in [3J [23] between the hydrodynamic large deviation functional 
and the one associated to the stationary measure: a density deviation in the stationary measure can 
be interpreted as the optimal space/time cost to produce this deviation starting from the steady 
state. Thus, one could try to see whether the optimal trajectory to generate a given steady state 
fluctuation in the WASEP converges to the optimal trajectory in the TASEP (after all, if the large 
v limit of the WASEP describes the TASEP steady state at the level of the density large deviation 
functional [20], one might expect the large drift limit to hold at the level of the whole trajectories). 

Our calculation of the large deviations of the integrated current relies on the hypothesis that 
the optimal density profiles are time-independent and that the measure conditioned on a given 
current is in local equilibrium. This led in particular to (|5.2|) which we do not expect to hold in 
the TASEP: as for the ring geometry, the measure conditioned on such an increase of the current 
should be significantly different from a Bernoulli measure. 

Even if we think that the current large deviations in the TASEP are correctly given by (|3.9|) . 
(jXTU)) . (j3~Tfij) (i.e. the limit v -> oo of the WASEP), we believe that the distribution of the 
fluctuations of the TASEP are not of diffusive nature and cannot be understood by taking simply 
the large v limit (for example in the WASEP the fluctuations of density are Gaussian |12| whereas 
in the TASEP they are not in the case of open boundary conditions |13|). 

Finally, it would be interesting to extend our approach to particle systems with more general jump 
rates. Hydrodynamic equations were derived for a broad class of microscopic models [2111261 l2"5] . In 
this case, there exists a density dependent diffusion coefficient D(p) and a conductivity coefficient 
a{p) such that (|2.1U|) becomes 

V(x,t) £»x [0,3], d t p = V x (D(p)V x p)-vV x [a(p)) . (6.1) 

The counterpart of the large deviation functional (|2.14| can be expressed in terms of the coefficients 
D and a [HI 13 El 

q-u f ■ \ f T M f 1 , Q/(M) -v<r(p{x,t)) +D(p(x,t))V x p(x,t)) 2 

%T] (J, p) = j Q dt ^ dx . (6.2) 

One can address the same type of questions as in this paper and minimize this functional under 
a current deviation constraint. The limiting expression (obtained for large drift u) might provide 
the extension of the Jensen- Varadhan functional to more general asymmetric dynamics. 

7. Appendix 

In Section 13.21 a constant flux was imposed and we proved the convergence of the large deviation 
functionals for v large. We are now going to investigate the effect of a time dependent current 
constraint. More precisely, we are looking for the optimal way of producing a current with different 
values on both sides of a singularity evolving at velocity v 

[ vcr(a), if x <C vt, 
I zat(o), if x 3> vt . 

For simplicity, we consider the hydrodynamic evolution in M instead of the finite system B = [0, 1]. 
The large deviation functional introduced in (|2.14|) can be defined in M. We would like to minimize 
the functional IY Q T , under the current constraint (|7.1|) and find the optimal profile p and the flux 
j (both depending on v) for v large. By construction one has 

d t p = -Vj . 
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Furthermore, we restrict the class of density profiles to the traveling waves so that 

p(x, t) = p(x — uvt), and vvp' = j . 

Integrating the previous equation, we deduce that j = v(yp + c) and the constraint (|7.1I) imposes 
that 

a (a) — a(b) aa(b) — ba(a) 

v = ——^ and c = — — — — = -va + a (a) = -vb + a(b) . 

a — b a — b 

We remark that this imposes the same condition on the velocity as in (|4.5j> . 
The large deviation functional minimized over the traveling waves is given by 

Gi,,a, b) = Ita inf/l fit f dx il^l±Ml\ (7.2, 

(v{c + V P )-V(T(p) + yf \ 

2 a(p) J 
where we used the explicit form of j in terms of p. 

We are going to recover the convergence to the Jensen- Varadhan functional (|4.8|) . The optimal 
cost is given by the profile p which minimizes (|7.2|) . We decompose G into two terms 

2J» <r(p) MP) Jm MP) 1 ' 

The functional G2 does not depend on the profile p 

1 b 




clog—^ p - v log(l - p) 

1-/9 



2j a P \a(p) a(p)J 2 

Using the explicit formula for c, we get 

G2 = \ [g(b) - g(a) - v{h(b) - h(a))] = ~F(a, b) . 

It remains to determine the optimal profile p by minimizing G\. Following the same computation 
as in we see that there is a constant K such that 

2 



\A =v 2 {c + vp-a{p)f + 2Ko(p). 



The constraint (|7.1jl at infinity implies that K = 0. Thus we get 

Gl = U dx ( c + ^-^)V =G2 
2 Jr cr(p) 

Combining the previous results, we get 

uF(a, b), if a > b 

0, if a < b ' 



G{vq) 



For v large (and in fact for all u), G{yq)jv is equal to the Jensen- Varadhan functional l|4.7|l . (|4.8|) . 
As expected the shocks (a < b) behave differently from the anti-shocks (a > b). The optimal 
profile has a sharp slope in an interval of size 1/v centered around vt and the cost to maintain the 
anti-shock is essentially located in this region. This confirms the Jensen- Varadhan theory which 
asserts that the contribution of each anti-shocks decouple (|4.1()j) . For this reason the previous 
approximation procedure should be also valid for a more general deviation of the current with 
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several discontinuities: in the limit v — > oo, each anti-shock should contribute independently to the 
large deviation functional. 

In jU], it was argued that some current deviations for the WASEP on a ring are achieved by 
a traveling wave. In the large drift limit, the exponential cost of these traveling waves was also 
related to the Jensen- Varadhan functional. 

Acknowledgments. We are very grateful to C. Bahadoran for explaining us the Jensen- Varadhan 
theory and for very useful discussions. 
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